Cooper-pair coherence in a superfluid Fermi-gas of atoms 
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f^**) , We study the coherence properties of a trapped two-component gas of fermionic atoms below the 

BCS critical temperature. We propose an optical method to investigate the Cooper-pair coherence 
CN| ' across different regions of the superfluid. Near-resonant laser light is used to induce transitions 

between the two coupled hyperfine states. The beam is split so that it probes two spatially separate 
• regions of the gas. Absorption of the light in this interferometric scheme depends on the Cooper-pair 

| coherence between the two regions. 

CO 

J* ; I. INTRODUCTION 

The coherence properties of clouds of trapped bosonic atoms have been already the focus of several interesting 
experiments. It has been convincingly demonstrated that the single-component Bose-Einstein condensate is coherent 
, on the time-scale of most experiments Q , while for two-state fluids the intra- and inter- component coherence survives 
on remarkably long time-scales despite external manipulations which displace spatially the centers of mass of the two 
Q , fluids with respect to each other [Q]. 

In the case of trapped two-state Fermi gases with negative scattering length, the existence of a BCS transition has 
been predicted ||. By now, intensive experimental work on trapping and cooling fermionic atoms Q with lower and 
lower fractions of Fermi temperature being reported and the success in trapping mixtures of bosons and fermions as 
O \ well as two states of the same fermionic atom allow most of the researchers to be optimistic that the BCS critical 

temperature will soon be reached. 
7-H ', The main candidates for atoms to undergo the BCS transition are 40 K and 6 Li Q. For 6 Li, the states |5) and |6) as 
defined in ^ have an anomalously large negative s-wave scattering length - thus a relatively high critical temperature. 
They have also a small enough decay rate in the region of mechanical stability, which gives the system a long enough 
lifetime to be experimentally interesting ||. It is also assumed that the densities of the two components are close 
enough so that below the critical temperature one has the usual s-wave BCS coupling between the two states. We 
will denote the two states involved in pairing by | J.) and | |). The BCS transition is connected with the appearance 
of an order parameter Q A(r) = — 4ti ^ a (ip^ (rjip^r)}. The ultraviolet divergences associated with using a contact 
interaction pseudopotential can be dealt, in case of metallic superconductors by introducing a cut-off at the Debye 
energy, and in case of trapped atoms by other renormalization schemes, such as pseudopotentials regularized with the 
operator <9 r [rx] The properties of the order parameter of cold fermionic gases in the superfluid phase have been 
I , studied in the regime close to the critical temperature, where it satisfies a Ginzburg-Landau equation |Q. 

In this paper we study the spatial coherence and pairing properties of the BCS-paired two-state Fermi gas at 
low temperatures, when the GL treatment is not applicable but instead one needs to use the Bogoliubov-deGennes 
formalism We propose a method for checking experimentally the coherence of Cooper pairs and measuring the 
Cooper pair size. The method is based on breaking some of the Cooper pairs by driving transitions between the 
states | t) and | J.) with nearly-resonant light. The laser beam is split and the two resulting beams are focused on two 
spatially separated regions of the condensate, after which they are merged together like in a typical intcrferometry 
experiment. If both the laser and the condensate coherence are preserved between the two regions, interference should 
have an effect on the absorption of the light. We show that this is indeed the case and the existence of interference 
contributions for a given beam separation shows that the gas has coherence on that length scale. 

Interaction of the superfluid gas with laser light has been considered earlier as a method to probe the BCS transtition 
of the gas P,pO|. In these proposals the paired states were coupled to some unpaired states of the atoms, and in the 
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proposals M the light was mainly thought of being off-resonant. Here we consider laser- induced transitions directly 
between the two states that are paired. The light is nearly resonant and supposed to be absorbed in the pair breaking 
process. Moreover, the interferometric configuration in our proposal allows to probe coherences, not only the existence 
of the superconducting gap. 

We first describe the proposed method in section 2, then derive the absorption rate for a general laser geometry in 
section 3. Absorption rate and the effect of coherences in the interferometric setup are considered in section 4. The 
appendix discusses coherence length scales in the BSC Green's functions. 

II. THE IDEA 

The basic idea of coupling the two paired states by a laser is depicted in Fig.l. The two states | f) and | J.) are 
typically hyperfine ground states thus one would actually use a Raman transition to couple them. In the following we 
will, however, consider just laser light with one Rabi frequency and detuning - in case of a Raman transition these 
can be understood as effective quantities. Likewise, Fig.l also shows only one laser driving the transition. 

One may make an educated guess that the laser has to have a finite detuning S to drive the transition. This is 
because one has to break the Cooper pair, i.e. provide extra energy at least of the amount of twice the gap energy 
in order to transfer and atom from, say, state | !) to state | f) which makes both of the pairing partners to have the 
same internal state and thus become excitations in the superfluid. Our calculations show that this extra amount of 
energy is indeed required. 

The scheme for probing spatial coherence is shown in Fig. 2. The laser beam is split and focused into two regions 
with a separation z. After this they have to be recombined, or detected nonselectively - this prevents from gaining 
information about from which of the beams light was absorbed. This lack of information is essential for any interference 
phenomena to occur. Below we calculate, first for a general laser configuration, the transfer rate of atoms from one 
state to another - which directly corresponds to the rate of absorption. Then we demonstrate the idea by considering 
a specific interferometric configuration. 

III. THE ABSORBTION RATE - GENERAL CASE 

In the rotating wave approximation the interaction of the laser light with the matter fields can be described by 
a time-independent Hamiltonian in which the detuning d plays the role of an externally imposed difference in the 

chemical potential of the two states. The total Hamiltonian becomes then H = H + Ht, where 

H = H B cs + U + ~) J dr^\ (r)^ (r) + L - J dr^\ (r)^ (r) . (1) 

Here /i is the chemical potential of the Fermi gas before the laser was turned on, the same for both of the components 
in order to allow standard BCS pairing. The Hamiltonian Hbcs is the BCS-approximation of the matter Hamiltonian 
with the chemical potential included . The transfer Hamiltonian is given by 



H T = J dm(f)i>] (ftfl (r) + n* (r)V] (f)V>t (f), (2) 

with f2(r) characterizing the local strength of the matter-field interaction. We consider that at the time — oo the system 
was in its BCS ground state, ^>bcs- Then the laser field has been turned on. The intensity of the electromagnetic 
field is small enough for the transfer Hamiltonian Ht to be just a perturbation on Hbcs- 

The main observable of interest is the rate of transferred atoms from, say, state J, to state f- This also directly 
corresponds to the absorption of the light. It is defined 

/t = Ft J rf «)|^(r> T (rl|*(i)) (3) 
and can be further evaluated with the help of the Schrodingcr equation ih-^\^!(t)) = H\^(t)) as 

J T = i [ dr(*(t)\n*(r)rp\(r)^(r) - fi(f^}(f)^(f)l*(t)>- ( 4 ) 
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In the following we call If the current in analogy to metallic superconductors where the flux of electrons out of the 
superconductor constitutes the electrical current. 

We introduce an interaction representation with respect to H and use linear response theory with respect to Ht 
jllj . Validity of the linear response theory requires that the laser intensity is small and the transfer of atoms can be 
treated as a perturbation. This implies also that the amount of atoms transferred from one state to another is not 
so big that it would severely imbalance the chemical potentials of the two species of atoms and thus break the BCS 
state. In order to simplify notation we denote 



6{t) = e 

where O is any operator. The current becomes 



iHBCStQg — iHBCSt 



(5) 



h = df 



BCS\ 



j(r,t),H T (t') |* 



BCS)- 



(6) 



Here for example 

H T (t) = J dr^^rje-^^jl^t^^^ + ni^^t^^t), 
j(f,t) = n*(f)e~ t5t ^l(r,t)^(f,t) - n(fOe i5 ^}(f,t)^(f,t) 

are the transfer Hamiltonian and the current operator in the representation (|^). We introduce the operator A(f, t) 
VL*(f)^yAr, t)^(f, t); using this the current of atoms to the f state becomes 

It = 



dt > {e iSit'-t) J df j df ^^ BCS \ [i(r,t),it(rV)J \^bcs) 
e «(*-f) f df [ df(^ BCS \ \At(r,t),A(ft,t')\ \^ B cs}}- 



This expression can be calculated by introducing a Matsubara correlation function 



X{iu>) 



dre iuJT / df / df(T T A(f,T)A^(f,0)) 



(7) 
(8) 

(9) 



which gives the rate of transfer by If = — 21 m [X ret {— 5, fj\ with the retarded correlation function defined by 

X ret (uj) luJ ^~ +l£ X{iu>). In the spirit of Wick's theorem, we split the Matsubara function into a part which con- 
tains only the superconductor Green's function Q and a part which contains only the anomalous superconductor 
Green's funtion T: X = Xjr + Xg, 



Xg(iu) = J dre luJT J df J 



df / dfn*{f)n{f)g{f,0;f, T )g(f,T;f,o) 



X?{iuj) 



dre 1 



df I df*n%f)n(f')J r \f,T-,7,Q)J r {f f ,T;f',0). 



(10) 
(11) 



After summation over Matsubara frequencies, and choosing zero temperature and 5 > 0, we get If = If .g + If ^ with 



I T , g = -27r^ / dfU(f)v n (f)u m (f) S(e n + e m - S) 

n,m 

I T ^ = 27T^ j dfdfn*(f)n(f)u* n (f)u m (f)v; n (f)v n (f)S(e n + e m ^S). 
Here the triplet (u ni v n );e n is a solution of the (nonuniform) Bogoliubov-de Gennes equations 



-e n v n (f) 



^V 2 + V(f) + U(f) - fi 
2m 



— V 2 + V(r) + U(f) - n 
2m 



u n (f) + A(f)v n (f) 
v n (f) - A*(f)u n {f). 



(12) 
(13) 

(14) 
(15) 
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Here V(r) is the external potential and U(r) = (^] (r)^ (r)) = ^^(tt>l (r)4>l(r)) is the Hartree field. Given 

the potential V{r) one can solve the Bogoliubov-deGennes equations and insert the results into (|l^)-([l3|). 

^From the equations (|l2|)-([l3|) on can immediately see that given a non-trivial spatial dependence of the laser profile 
f2(r), the BCS coherence, that is, the spatial coherence of v n (r)u m {r) in different length scales will have an influence 
on the current. In the following section we demonstrate this by a simple example. 



IV. THE ABSORPTION RATE - INTERFEROMETRIC CASE 



We consider the setup of Fig. 2, where the laser beam is split and focused in the z-direction but homogeneous in 
the two other directions, that is Sl{r) = £1 [S(z — z{) + S(z — z-z)]. This gives stronger ^-dependence than focusing into 
two points, and is probably also experimentally simpler. 

In order to derive instructive results in closed form we consider the case of a uniform gas. This corresponds to 
working in a region of space much smaller than the oscillator length of the trap, or having a shallow trap where 
confinement effects are small. Solutions of the Bogoliubov-deGennes equations are given in this case by 

u k (r)=u k e^ ; u\ = \ (\ + , (16) 

Vk{f)=Vke ikr . ,| = lfl-|j, (17) 

^ = V^ 2 + A2 ' ^ISST""- (18) 

Cylindrical coordinates and momenta are used and the sums are transformed into integrals. The currents (|l2])-(|T^) 
become 







\n 



h.Q = —t / uT X t 1 + cos ( fc * + 5(E q + E k - S) (19) 



h,r = -y I u k v ku q v q [1 + cos(fc z + q z )z] 5{E q + E k - 8). (20) 



Here k = p 2 + k 2 , q = p 2 + q 2 , z = z\ — z-i , and the integral symbol means J = J Q dpp J Q dk z J Q dq z . One can 
investigate how absence of coherence affects the results by attaching to u k {f) and v k (f) phase factors which describe 
random space- and time-dependent fluctuations and by averaging over them. This makes the cosine-dependent term 
in I^g as well as the whole current to disappear. Thus one can vary z and see whether the current varies as well 
— if the current becomes constant for large enough z one knows the oscillating terms are absent, that is, coherence 
is not preserved in that length scale. 

The currents I^g and I^jr are presented in Fig. 3 and 4. In the numerics the 5-functions in (fL9|)-(po|) have been 
replaced by Lorentzians of width 0.01 to describe the finite laser linewidth. It is interesting to compare the variation 
of the currents as a function of z to the typical length scale of Cooper pairing, the Cooper pair size xqp- The pair 
size can be estimated simply by the momentum-coordinate uncertainty relation, the uncertainty in momentum being 
determined by the corresponding uncertainty in energy (the gap). This results to xqp — jv/?. In the case of Figs. 3 
and 4, xcp ~ 7. One observes oscillations occuring with approximately the period Az ~ 2 which is of the same order 
of magnitude as xcp- For z > xcp the current gradually becomes constant. The total current 1^ g + is shown in 
Fig. 5. Exactly at S = 2A the currents I^g and Jf y cancel each other (this would be different in a non-homogeneous 
case) but for larger values of detuning the current is non-zero. The oscillations are visible. 

Another length scale involved is the Fermi length xf given by xf = l/\/2mp. The relationship to the Cooper 
pair size is typically xcp ^> xf, in case of Figs. 3 and 4 ij? ~ 0.7 = O.l^cp. In physical units, xcp ~ 0.3/jto 
and xf ~ 0.03/im in a system with particle number N ~ 10 7 and trap frequency — 2ir x 150Hz (and the gap 
A = O.lEp)- If N ~ 10 5 , the length scales are xcp ~ 2pm and xp ~ 0.2pm. Thus for realistic experiments, the 
Fermi length will always be below the diffraction limit of light and cannot be resolved. This is the case for xqp too 
for very high particle numbers and/or trap frequencies. But for moderate densities - which is expected to be the 
situation in the first experiments - the Cooper pair size may be of the order of several micrometers which can be 
resolved by light in principle. 

Note that the detuning S has bigger than twice the gap to produce any current at all (there is some current below 2A 
because the 5- functions in (|l9|)-(p0|) were replaced by Lorenzians). Thus having or not having BCS pairing (nonzero 



4 



or zero gap) makes a significant difference in the result, which means that the method can be used also to detect 
the onset of the phase transition, c.f. Jo| , [Tc| , |T2| . |T3| | . For certain types of perturbations the Cooper-paired atomic Fermi 
gas has also below-gap excitations. A density perturbation term of the form U(r, t)[^j(r)^-f (r) + $t (r)-0j (f)\ leads 
to appearance of a Bogoliubov- Anderson phonon |l2). In our case the laser detuning gives a similar term in the 
Hamiltonian, S/2[^(r)'tpi(r) — ^{r)^ [{r)] (assuming that the laser is turned on and confined to a spatial region, one 
may think of 5 having space and time dependence). The minus sign, however, leads to cancellation of certain terms 
in the response calculation and to the absence of the usual Bogoliubov- Anderson phonon. This will be discussed in 
detail in another publication |fl4|. 



V. CONCLUSIONS 



For standard metallic superconductors the coherence properties are usually investigated by creating interfaces with 
normal metals or with insulators. Optical manipulation of cold alkali gases offers the unique opportunity of creating 
such interfaces at any point in space and of controlling at will the transfer Hamiltonian across the interface. We 
described a method for testing BCS coherence in a system of atomic fcrmions cooled below the critical temperature 
by driving laser-induced transitions between the two paired hyperfine states in different regions of the space. We 
presented numerical and analytical evidence for the feasibility of the procedure. We show that the photon absorbtion 
rate changes when the distance between the transition regions is of the order of magnitude of the BCS correlation 
length. Absence of coherence would make the rate independent of the distance. The method can serve as well for 
the detection of the onset of BCS transition, because the absorption peak is shifted to the detuning of twice the gap 
energy. 
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VII. APPENDIX 



The currents I^g and ij jr are convolutions of pairs of correlation functions, thus their direct connection with xcp 
and Xf is not transparent. In this appendix we discuss the dependence of the BCS correlation funtions on these length 
scales. It gives insight in the physics of Cooper-pair coherence and to the more complicated length scale dependence 
of I^g and ij,jr. 

The Green's funtions for the superconductor are defined as 



T{r,T;r,T) = -{%■ ^ T (r, r)ijj L (f, t') ) 
T\r,T-7,T') = -{% [$(?,T)$(?y)]) 
g(r,T;r , ,T , ) = -{T r ?/> T (f,r)Vi|(r,T') ). 

In the uniform case and at zero temperature we get for the zero-time correlation functions 



(21) 
(22) 
(23) 



m 
47r 2 r 

m 
An 2 r 



dt; sin y/2m(n + £)r 
<i£ sin \/2m(/i + £)r 



1 - 



A 



(24) 
(25) 



where r — \r — f*|. Information about any possible coherence length scales like xcp or xf = l/\/2m/i is hidden in 
these equations. Note that for a normal system (zero gap) the Fermi correlation has quite an explicite dependence on 
Xf- 
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G(f, 0; f*, 0) A ==° -t— l±E. sin J- - ±E. cos — . (26) 



a=o 1 [ xi, . r xt r 
— — -3- sin J- cos — 

2it z xf \ r- 3 xf t* xf 

The anomalous correlation function T is non-vanishing only when the gap is finite, reflecting the instability of the 
Fermi sea and the appearance of fermionic pairing. For T one can actually derive an approximate formula which shows 
the dependence on xqp or xf- For the usual situation A <C M, and one can expand y/2m(jI+£) 



r 



which gives 



1 r g 

Kjf 1 2 icp 



^0;^0)_sin^-jX ^). (27) 



Here Kq is the zero order modified Bessel function which for a real argument a admits the integral representation 

cos(ax)dx 



K (a) 



VT + r 



2 



For example, using the series expansion of Kq we find that at r small compared with xcp the anomalous correlation 
function decreases with r as ha.(r/xcp)/r, 

toA sin — ( r \ 

F(f,0;?,0) = = —— in- + 0.577215). (28) 

2-K l r \ ixcp J 

In conclusion, the spatial behavior of the normal and anomalous Green's functions is governed by two parameters, 
xcp and xp- However, since Xf is typically much smaller than the width of the laser beams we expect a dependence 
only on xcp in the result of a measurement. 
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FIG. 1. Laser coupling between the two paired states. Atoms in two different hyperfine states | f) and | J.) have a negative 
s-wave scattering length and form a Cooper pair associated with the appearance of the superconductor order parameter (gap) 
A. The laser couples the states | f) and | J,), the frequency of the laser is detuned from the atomic transition frequency by the 
amount 8. In order the transfer of the atom from one state to other to happen the Cooper-pair has to be broken and the atoms 
become excitations in the superfluid. Note that in a real experiment a Raman scheme would be used instead of a single laser 
beam. 

FIG. 2. Schematic picture of the interferometric setup for probing coherences. The initial laser beam is split coherently and 
focused on two spatial regions of the gas with the separation z. After the beams have passed the gas they are recombined again 
and the amount of absorption is measured. 

FIG. 3. The current I^ t jr as a function of the separation z and the detuning S. The current oscillates with z, with the period 
of ~ 2 ~ xcp/3. The current is peaked around S — 0.2 = 2A. Here the Fermi energy /i = 1 and the gap A = 0.1. All the 
variables are dimensionless. For instance for the particle number N = 10 5 and the trap frequency Q — 2-k x 150Hz the scales 
would be [z] ~ 0.3/im and [5] ~ 75kHz. 



FIG. 4. The current —I\g. Oscillations similar to those in I^j: are observed. The parameters are like in Fig. 3. 



FIG. 5. The total current / = \I^g + Oscillations are observed as in Figs. 3 and 4. 
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